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Abstract. In this article we review and discuss a mechanism for coupling between 
electric polarization and magnetization that can ultimately lead to electric-field 
switchable magnetization. The basic idea is that a ferroelectric distortion in an 
C$ ' antiferromagnetic material can "switch on" the Dzyaloshinskii-Moriya interaction 

^ which leads to a canting of the antiferromagnetic sublattice magnetizations, and 

| thus to a net magnetization. This magnetization M is coupled to the polarization 

£H • P via a trilinear free energy contribution of the form P ■ (M x L), where L is the 

antiferromagnetic order parameter. In particular, we discuss why such an invariant 
is present in R3c FeTi03 but not in the isostructural multiferroic BiFe03. Finally, 
we construct symmetry groups that in general allow for this kind of ferroelectrically- 
induced weak ferromagnctism. 
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1. Introduction 



> 
ON 
00 

in 
o 

00 \ One of the key challenges in the field of multiferroics is to design and/or discover 
materials that exhibit a strong coupling between magnetic and ferroelectric order 
parameters. Such materials are of fundamental interest as they provide a novel platform 
to study how microscopic degrees of freedom, such as spin and lattice, interact to produce 
macroscopic phenomena. These strongly coupled multiferroics are also anticipated 
to find application in future generations of novel devices in which magnetization can 
be controlled via an electric-field and/or electric polarization can be controlled via a 
magnetic field. 

Two different ways of controlling the state of a magneto-electric system are possible: 
phase control [H [2j [3] or domain control [HE]. In the first case an external field is 
used to trigger a phase transition between two fundamentally different phases. By 
tuning to the vicinity of a phase transition where two such phases compete |6j [7], 
e.g., an antiferromagnetic-paraelectric phase and a ferromagnetic- ferroelectric phase, a 
large magneto-electric response can be produced [8] even in systems where the intrinsic 
coupling between magnetic and ferroelectric order parameters is not very strong [9] . In 
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the case of domain control, the external field triggers a transition between two equivalent, 
but macroscopically distinguishable domain states, i.e., different realizations of the same 
phase. Here, the magneto-electric coupling has to be large enough to overcome the 
energy barrier for domain switching, which in general depends on both the initial and 
the final domain. Once the field is removed, the system is stable in the new domain 
state. 

Several materials have been identified that realize one of the two scenarios described 
above. In particular, much recent research has been focused on systems where magnetic 
order itself breaks spatial inversion symmetry and electric polarization therefore appears 
as a secondary order parameter [HEJ ED]- Several different microscopic models have been 
proposed that lead to such "magnetically-induced ferroelectricity." In some cases this 
effect is caused by the presence of spin-orbit interaction [HI [12], whereas in other cases 
spin-orbit coupling is not required [HH E] • See other contributions in this Focus Issue 
for a more detailed discussion of this interesting topic. 

In the present article we discuss a somewhat different, but in a certain sense 
complementary, possibility to realize coupling between magnetic and ferroelectric 
order parameters, which can then be used to achieve domain control of the 
corresponding multiferroic system. The basic idea is that a ferroelectric distortion 
in an antiferromagnetically ordered material can cause a small magnetization due 
to the Dzyaloshinskii-Moriya interaction. In this case of "ferroelectrically-induced 
ferromagnetism" [15] the polar distortion gives rise to both the electric polarization and 
the magnetization, hence the two quantities are inherently coupled. Here, in contrast 
to the case of magnetically-induced ferroelectricity mentioned above, the magnetization 
is the secondary order parameter that is coupled to the primary order parameter, the 
electric polarization. 

The concepts and ideas reviewed and discussed in this article are based mostly 
on a series of publications by the present authors, and recently have been used to 
identify a specific class of materials that are predicted to exhibit this effect [T6 | [T71 [18] . 
The concept of ferroelectrically-induced ferromagnetism was first suggested by Fox and 
Scott [15] based on macroscopic symmetry properties for the magneto-electric fluoride 
BaMnF4. Here we discuss a specific microscopic mechanism leading to such macroscopic 
behaviour and analyze the corresponding symmetry requirements. 

In the following we first review the basic idea behind the proposed mechanism, then 
summarize our previous work on this topic. We choose to illustrate the general concept 
by discussing one structure in detail, the ten atom rhombohedrally distorted ABO3 
perovskite structure, paraelectric space group R3c, and the corresponding ferroelectric 
subgroup R3c, although the established principles are easily generalizable. Specifically 
we show why the proposed effect is present in magnetic A site R3c perovskites such 
as FeTiC>3, but not in R3c BiFeOs, where the magnetic cations are situated on the 
perovskite B sites. Finally, we present a rather general discussion of the various 
symmetry aspects that have to be taken into account in order for a material to exhibit 
the desired behaviour. Ultimately our goal is to highlight the unique and powerful 
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approach of combining effective microscopic models with symmetry arguments to guide 
first principles calculations in the discovery of new phenomena and the design of their 
material realizations. 

2. Weak ferromagnetism and electric polarization 

It was realized by Dzyaloshinskii in 1957 that the appearance of "weak" ferromagnetism 
in some antiferromagnetic materials such as e.g., Fe2C>3, and its absence in the 
isostructural system Cr 2 03 can be explained solely on grounds of symmetry [19]. The 
symmetry of a magnetically ordered material depends on the underlying crystallographic 
structure, the orientation of the magnetic moments relative to each other, and on the 
orientation of the individual magnetic moments with respect to the crystallographic 
axes [201 l2Tj . Dzyaloshinskii showed that an invariant in the free energy expansion of 
the form 

E DML — D • (M X L), (1) 

— * — * — * 

where D is a materials-specific vector coefficient, M is the magnetization, and L 
is an antiferromagnetic order parameter, results in the secondary order parameter 
M appearing at the antiferromagnetic ordering temperature. In other words, if the 
symmetry of the purely antiferromagnetic state is such that the appearance of a small 
magnetization does not lead to a further symmetry lowering, then any microscopic 
mechanism which favours a nonzero magnetization, even if it is rather weak, will lead 
to M ^ 0. 

It was later shown by Moriya that an invariant of the required form can result from 
an antisymmetric microscopic coupling between two localized magnetic moments Si and 
S f . 

ET = 4 • (3 x Sj) , (2) 

and that such an interaction arises from the interplay between superexchange and spin- 
orbit coupling [22]. Invariance of the interaction energy (F2J) with respect to exchanging 
Si and Sj requires that cfy = — dji. The energy ([2]) is minimized when the two magnetic 
moments form a 90° angle (or more accurately when d^, Si, and Sj form a left-handed 
system for \dij\ > 0), but due to the simultaneous presence of the generally much stronger 
Heisenberg-type interaction Ef, = JijSi ■ Sj, with = Jji, which favours either or 
180° angles, the presence of the Dzyaloshinskii-Moriya (DM) interaction usually only 
leads to a small canting between the interacting moments, i.e. a small deviation from 
an overall collinear magnetic configuration as illustrated in Fig. [TJ 

Symmetry puts strict constraints on the DM vector dij. In particular, d^ is 
identically zero if the midpoint between the magnetic moments is an inversion center. 
Turning this around, it means that in magnetic ferroelectrics, which do not exhibit 
any inversion centers, the DM interaction can be expected to be a rather common 
phenomenon (even though there can be other symmetry restrictions which prohibit 
the DM interaction even in the absence of any inversion centers). This raises the 
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Figure 1. The presence of the Dzyaloshinskii-Moriya interaction (D ^ 0) leads to a 
slight canting of the magnetic moments and a resulting net magnetization M ^ in 
an otherwise collinear antiferromagnet. 



question of possible cross-correlations between electric polarization and DM interaction 
in ferroelectric magnets. 

Many ferroelectric systems of interest today can be characterized by a small polar 
structural distortion away from a centro-symmetric reference structure. This reference 
structure is usually identical with the crystal structure of the paraelectric phase but 
need not be [23]. The distortion, which leads to an electric dipole moment, can be 
reversed by applying an appropriate electric field (domain-switching). This leads to 
the following scenario [HU [18]: If in the paraelectric reference structure the midpoint 
between two neighbouring magnetic ions is an inversion center, which is destroyed by 
the ferroelectric distortion, then the ferroelectric distortion can "switch on" the DM 
interaction between the two ionsjj] In this case the materials specific parameter, D, of 
Dzyaloshinskii's invariant, Eq. (TO , can be identified with the electric polarization P and 
the corresponding invariant in the free energy expansion is: 

E PML ~ P ■ (M x L) . (3) 



From Eq. ([3]) it can be seen that if it is possible to reverse the direction of P (using an 
electric field) without changing the orientation of L, then the magnetization will reverse 
too, in order to minimize the total free energy of the system. Reversal of L can be 
prevented by sufficiently large magnetic anisotropy, and thus an invariant of type ([3]) 
opens up the possibility for electric field-induced magnetization switching. 

In the preceeding paragraph we have outlined the general scenario that can lead 
to electric field-switchable weak ferromagnetism. In order to find specific example 
materials that exhibit this effect it is important to point out that even for cases where the 
ferroelectric distortion destroys the inversion centers between adjacent magnetic sites, 
and thus provides a necessary requirement for nonzero DM interaction, there can be 
other symmetry operations that result in dij = or prevent the system from exhibiting 
a macroscopic magnetizationj§] In the following we quickly summarize our previous 
work along these lines and then discuss in detail why the desired coupling is present 
in R3c FeTiC>3 but absent in BiFeOs, even though these two systems have the same 
crystallographic space group symmetry. 

| This has been called the "structural-chemical criterion" in Ref. [18] . 
§ This has been termed the "magnetic criterion" in Ref. [18] . 
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3. First principles studies of BiFe0 3 , BaNiF 4 , and FeTi0 3 

The possibility of an electric-field switchable DM interaction was first discussed in 
the context of BiFeO"3 [16]. BiFeC>3 is an antiferromagnetic ferroelectric with a Neel 
temperature of ~ 643 K and a ferroelectric Curie temperature of ~ 1103 K [24l [25] . It is 
thus a very rare example of a multiferroic with both magnetic and ferroelectric ordering 
temperatures above room temperature, and it is probably the most-studied multiferroic 
to date. The primary magnetic order in BiFeC>3 is "G-type" antiferromagnetism |26j, 
i.e., "checkerboard" -like, but in addition it has been reported that bulk single crystals 
exhibit a superimposed cycloidal spiral magnetic ordering with a large period of ~ 620 
A [27]. This spiral ordering seems to be absent in thin film samples [28] where instead 
a small magnetization has been found [29l [30] [28] . Early reports of magnetizations up 
to l/xe/Fe [29] could not be confirmed in other samples [30] 28J and therefore seem to 
be caused by extrinsic effects. 

In [16] it was shown by first-principles calculations that spatially homogeneous 
BiFeO"3 (without the spiral spin structure) exhibits weak ferromagnetism as a result 
of the DM interaction, and that the resulting magnetization is about 0.1 /ie/Fe. 
Furthermore, it has been shown in [31] that the magnetization is only weakly dependent 
on epitaxial strain, another indication that the very large magnetization of about 1 
/xe/Fe found in the original thin film samples |29j is most likely due to extrinsic effects. 
The calculated magnetization agrees well with more recent experimental thin film data 
[301 [28] . We point out that a meaningful comparison between measured and calculated 
magnetization requires that an antiferromagnetic mono-domain state has been achieved 
in the experiment. 

It was also shown in [16] from explicit first-principles calculations that the sign of 
the DM vector D in BiFeO"3 is independent of the polar distortion, but that it is instead 
determined by a rotational (non-polar) distortion of the oxygen octahedra network 
present in BiFeCV As we review in the next section, this lack of coupling between the 
DM vector and the polarization in BiFeO"3 is a question of symmetry [16] [18] and not 
due to a "weak coupling" as some authors have suggested [321 [33] . Despite this absence 
of coupling between P and M in BiFeO"3, it was nevertheless realized that the DM 
vector can indeed couple linearly to a structural distortion [16] . However, the relevant 
structural distortion in the case of BiFeC>3 (the octahedral rotations) is non-polar, and 
therefore it was concluded that electric-field switching of the weak magnetization in 
BiFeO"3 is unlikely as there is no obvious way to couple the electric-field to the non-polar 
distortion. Still, it was suggested that if materials can be found where the polarization 
and the DM vector are due to the same structural distortion, then electric field-induced 
reversal of the weak ferromagnetic moment is possible, and that there are no general 
symmetry arguments that prevent such an effect. A specific material that realized the 
predicted effect, however, remained elusive until the work of Ref. [IB] . 

A subsequent study of the antiferromagnetic ferroelectric BaNiF4 revealed that the 
DM vector can indeed be proportional to a polar ferroelectric distortion [IT]. However, 
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the overall symmetry of BaNiF 4 does not allow a macroscopic magnetization, and it was 
shown that even though the DM interaction leads to a local canting of neighbouring 
magnetic moments, all components cancel out when taking the sum over the whole unit 
cell such that M = 0. BaNiF 4 can therefore be classified as weak anti ferromagnet. It 
was pointed out that it should be possible to switch the secondary antiferromagnetic 
order parameter in BaNiF 4 using an electric field [TT] . 

Finally, it was recently shown in [18], that a series of R3c titanates, ATiC^ with 
magnetic cations A=Mn, Fe, Ni, do in fact combine all the symmetry properties 
necessary for ferroelectrically-induced weak ferromagnetism. These compounds can be 
synthesized at high pressure and remain metastable at ambient conditions [SUES]. It 
was confirmed in [18] by explicit first principles calculations that a weak magnetization 
exists in these materials, and that it is reversed when the polar distortion is reversed 
while keeping all other order parameters (apart from the magnetization) fixed. It can 
therefore be expected, that in these materials the magnetization can be reversed using 
an electric field. This represents the first specific example for the general mechanism of 
magneto-electric coupling outlined in Sec. [2] 

In [18J two criteria for the rational design of ferroelectrically-induced weak 
ferromagnetism were formulated: a "structural-chemical criterion" and a "magnetic 
criterion" . As already briefly remarked via footnotes in Sec. [2j the structural-chemical 
criterion implies that the midpoint between two magnetic sites is an inversion center 
in the paraelectric reference structure, whereas the magnetic criterion expresses the 
fact that there should be no other symmetry elements besides these inversion centers 
that prevent the system from exhibiting weak ferromagnetism. In the following we 
will present a detailed comparison between BiFeC>3 and FeTiC>3 (as representative for 
the titanate systems discussed in [18]), focusing on symmetry aspects, and discuss 
why in FeTiC>3 both criteria are fulfilled, i.e. the DM interaction is induced by the 
ferroelectric distortion, whereas in BiFeC>3 this is not the case. We will then develop 
guidelines to construct symmetry groups that generally allow ferroelectrically-induced 
weak ferromagnetism. 

4. Magnetic A-site versus Z?-site R3c distorted perovskites: symmetry 
aspects 

The R3c structure in which both BiFeC>3 and FeTiC>3 are found, can be regarded as a 
distorted version of the 5-atom cubic perovskite structure. The deviation from the ideal 
perovskite structure can be decomposed into two components (see Fig. [2]): 

I. antiferrodistortive counter-rotations of the oxygen octahedra around the [111] axis 
(leading to a unit cell doubling compared to the perovskite primitive unit cell), and 

II. polar displacements of all the ionic sublattices relative to each other parallel to 
[111]. 
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Figure 2. Schematic depiction of the structural distortions that lead from the perfect 
cubic perovskite structure to the R3c structure of BiFeOa/FeTiOs: I. counter-rotations 
of oxygen octahedra around [111] (red); II. ionic displacements along [111] (green) 

In FeTi0 3 the rotations (I) are so large that R3c FeTi0 3 is usually not considered 
to form a distorted perovskite structure but rather the "ferroelectric lithium niobate 
structure" jjj] Note that this LiNb03 polymorph of FeTiOs is structurally isomorphic 
to BiFeC>3 except that the positions of the Fe and Ti/Bi cations are exchanged, i.e.: 
FeTi0 3 -^BiFe0 3 implies Fe^Bi and Ti^Fe. 

Let us first consider BiFeO"3 and FeTiOs in the ideal perovskite structure, space 
group Pm3m. Here we neglect both the octahedral rotations (component I) and the 
ferroelectric displacements (component II) while assuming that the spins order in a G- 
type antiferromagnetic pattern. For both BiFeO"3 and FeTiOs this paraelectric reference 
structure has inversion centers at all midpoints between magnetic sites. Furthermore, 
these inversion centers will be destroyed by any ferroelectric distortion, and thus the 
structural criterion described in [18] is fulfilled for both systems. However, since the 
cubic perovskite structure contains only a single magnetic cation per unit cell, an 
additional symmetry operation exists in the antiferromagnetically ordered state that 
requires M = 0, i.e., the magnetic criterion is not fulfilled. The corresponding symmetry 
operation translates all ions by one unit cell and then inverts all magnetic moments 
through time inversion. In general, M = by symmetry, whenever the magnetic unit 
cell is a multiple of the chemical unit cell. Specifically for the case of BiFeO"3/FeTiO"3 
this implies that the unit cell doubling caused by the octahedral rotations (component 
I) is essential for obtaining weak ferromagnetism and to fulfill the magnetic criterion. 
Therefore, the octahedral rotations cannot be ignored in a proper symmetry analysis of 
BiFe0 3 or FeTi0 3 , and we conclude that cubic perovskite is not a suitable paraelectric 
reference structure for obtaining ferroelectrically-induced weak ferromagnetism. 

In contrast, if we only neglect the polar displacements (component II) while 
including the octahedral rotations (component I) in our paraelectric reference structure, 
we obtain the Calcite or "paraelectric lithium niobate" structure (space group R3c). 

| We point out that this distinction is not a question of symmetry but merely of structural connectivity. 
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Figure 3. Displacement vectors of the oxygen anions (red) corresponding to the 
octahedral rotations (I), relative to the ideal cubic perovskite structure. Due to 
the displacement of the oxygen anions located at the midpoints between adjacent B 
sites (green), the corresponding inversion centers are destroyed, whereas the inversion 
centers located at the midpoints between the A sites (blue) are conserved. 

This is the closest centrosymmetric reference structure for both BiFe03 and FeTiC>3 
and it is thus the proper starting point for a Landau free energy expansion describing 
the possible coupling between P, L, and M in these systems (see Sec. [5]). 

The crucial difference between BiFeC>3 and FeTiC>3 in the R3c structure is the 
different local site symmetry of the magnetic cations: in BiFeC>3 the magnetic cations 
(Fe 3+ ) are situated on Wyckoff positions 2b, corresponding to the B site of the underlying 
perovskite structure, whereas in FeTiC>3 the magnetic cations (Fe 2+ ) occupy Wyckoff 
positions 2a, corresponding to the perovskite A sites (see Fig. |2]). To see whether the 
structural criterion of [18] is fulfilled, we have to check whether the midpoints between 
the magnetic sites (i.e., the Fe sites) in each system are inversion centers or not. As 
illustrated in Fig. [3] this is the case for FeTiC>3 but not for BiFeC>3, i.e., the structural 
criterion for ferroelectrically-induced ferromagnetism is fulfilled in FeTiC>3 but not in 
BiFeC>3. Note that in the completely undistorted perovskite structure this criterion 
was fulfilled for both systems, i.e., inversion centers were located on the midpoints 
between the A sites as well as on the midpoints between the B sites. The latter are 
destroyed by the octahedral rotations (I). This is the reason why BiFeC>3 exhibits weak 
ferromagnetism already in the paraelectric R3c phase, i.e. weak ferromagnetism in 
BiFeC>3 is induced by the octahedral rotations and not by the polar distortion [16J. 

We point out that in some previous publications "i?3c" symmetry has been used 
as basis for a free energy energy expansion of BiFeC>3 (see e.g. [36j S3 E2])- "-R3c" 
is not standard space group notation^, and it is therefore not fully clear what the 

If Note that here the "bar" indicating space inversion is combined with the c-type glide plane and not 
with the threefold rotation as in R3c 
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corresponding symmetry operations are, but apparently in this case it is assumed that 
inversion centers are located between the magnetic sites. As stated above this is not true 
for BiFeOa in the proper paraelectric reference structure. It follows that the symmetry 
analysis in [361 [37J [32] in fact applies to FeTi0 3 instead of BiFe0 3 . 

It is now clear that only FeTi03, but not BiFeOs, fulfills the structural-chemical 
criterion. What remains to be discussed is whether the inversion symmetry related to the 
midpoints between the magnetic cations in FeTiC>3 is the only symmetry operation that 
requires M = 0, i.e., whether the magnetic criterion is satisfied. One candidate for such 
a symmetry operation was already mentioned earlier: the combination of a primitive 
translation of the paraelectric structure with time inversion. However, this can only 
be a symmetry operation if the magnetic unit cell is a multiple of the crystallographic 
unit cell. As discussed above, the presence of the octahedral rotations in R3c FeTiOs 
doubles the crystallographic unit cell compared to simple perovskite, and thus ensures 
that for the assumed G-type magnetic order the magnetic and crystallographic unit cells 
are identical. 

Apart from such "magnetic Bravais lattice translations" [21] the fulfillment of the 
magnetic criterion is mainly a question of how the individual magnetic moments are 
oriented relative to the crystal axes, which is determined by the magneto-crystalline 
anisotropy. As already shown by Dzyaloshinskii [19], an orientation of the magnetic 
moments perpendicular to the rhombohedral axis is required for weak ferromagnetism 
to occur in crystallographic R3c symmetry. First principles calculations show that this 
orientation is indeed favored by the magneto-crystalline anisotropy in FeTiOs [IE], and 
thus the magnetic criterion is fulfilled for this system. 

In summary, we have shown that in FeTiOs both the structural-chemical and the 
magnetic criteria for ferroelectrically-induced weak feromagnetism are fulfilled. Due to 
the different location of the magnetic cations in BiFeOs the structural criterion is not 
fulfilled in this system (and thus the magnetic criterion is not applicable). This shows 
that in FeTi03 the weak magnetization is linearly coupled to the electric polarization 
whereas this is not the case in BiFeOs, as was also verified by explicit first principles 
calculations in [16, 18J. It also becomes apparent that it is important to consider the 
full crystallographic symmetry to analyze possible coupling between P and M, including 
also nonpolar structural distortions such as the octahedral rotations in BiFeOs/FeTiOs. 

5. General symmetry considerations 

In the preceeding sections we described the general idea behind ferroelectrically-induced 
weak ferromagnetism, we reviewed results for some specific example materials, and 
presented a detailed comparison of the two R3c structure materials BiFeOs and FeTiOs. 
In Sec. [2] we mentioned that in order to achieve the desired coupling between the 
polarization and the magnetization, the free energy expansion of the high-symmetry 
phase has to contain a term of the form ([3]). In this final section we discuss some 
general symmetry aspects and we list magnetic point groups that are compatible with 
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the existence of such a coupling between P and M. 

According to the Landau theory of phase transitions, which describes continuous 
transitions from a high symmetry configuration into a configuration with lower 
symmetry, the free energy of the system can be expanded in powers of the various 
order parameters [381 139] . Each individual term in this expansion has to be invariant 
with respect to all the symmetry operations of the high symmetry phase. Here, we 
are concerned with a transition leading from a paramagnetic and nonpolar phase 
(where P = M = L = 0) into a lower-symmetry phase where all three order 
parameters are nonzero. This transition implicitly also defines an intermediate nonpolar 
antiferromagnetic phase where L ^ but P = M = 0. 

In order to achieve the desired coupling between P and M, the term P ■ yM x Lj 
has to be allowed in the free-energy expansion of the high symmetry phase. Since the 
transformation properties of the polar vector P and the axial vector M are well known, 
the question of whether such a trilinear coupling between P, M and L is allowed, 
is basically a question about the symmetry properties of the antiferromagnetic order 
parameter L. We point out that L cannot generally be classified as either axial or polar 
vector; its symmetry properties depend on the microscopic definition of L in terms of 
ionic magnetic moments and the underlying crystallographic symmetry. Ultimately, the 
design criteria of [18] aim at designing an antiferromagnetic order parameter with the 
required macroscopic symmetry to couple P and M. 

A first symmetry requirement for L has already been discussed in Sec. HJ the 
magnetic unit cell has to be identical to the crystallographic unit cell. Otherwise, 
a symmetry element exists which consists of time reversal combined with a lattice 
translation of the paramagnetic phase. Such a transformation leaves L and P invariant, 
but changes the sign of M and thus (j3J) is not an invariant of the corresponding high 
symmetry phase. This symmetry requirement for L is implicitly expressed in [18] as 
part of the magnetic criterion. 

The next symmetry operation to consider is space inversion. The presence of space 
inversion in the high symmetry group is not a necessary requirement , but it simplifies the 
following analysis considerably. We therefore restrict ourselves to cases where the high- 
symmetry paramagnetic, non-polar reference structure is centrosymmetric. Notice, since 
P is a polar vector that changes sign under space inversion and M is an axial vector that 
is invariant under space inversion, that the antiferromagnetic vector has to be odd under 
space inversion to allow an invariant of the form P ■ \M x Lj. It is apparent that such 
an antiferromagnetic order parameter arises for example if the inversion centers of the 
high symmetry structure are located between two antiferromagnetically coupled cations, 
but not if these inversion centers are located on the magnetic cation sites themselves. 
This symmetry requirement for L is thus related to the structural criterion discussed 
previously. 

We now construct symmetry groups that allow for an antiferromagnetic order 
parameter with the two symmetry requirements outlined in the two preceeding 
paragraphs. To simplify the presentation we only discuss magnetic point groups, not 
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the full magnetic space group symmetry. Note that the macroscopic point group 
corresponding to a particular microscopic space group is obtained by neglecting all 
translational parts of the corresponding space group operations [2T] . Therefore, the 

— # 

symmetry properties of L outlined above are equivalent to stating that the magnetic 
point group of the intermediate nonpolar antiferromagnetic phase can contain neither 
space inversion nor time reversal symmetries individually (i.e., neither 1 nor 1'), but 
only the combined operation of space inversion followed by time reversal (i.e., 1'). Note 
that these are also the symmetry requirements for the existence of a magnetic toroidal 
moment T (see [4"0| HI]), which hints at a close connection between the presence of a 
magnetic toroidal moment and the presence of the magneto-electric coupling expressed 
in ([3]) (the exact nature of this connection will be the subject of future investigations). 

In addition, we realize that a free-energy invariant of the form ([3]) with an 
antiferromagnetic order parameter that can be classified as a toroidal moment gives rise 
to an antisymmetric linear magneto-electric effect (see e.g. |41j). and we can therefore 
now list all the magnetic space groups of the desired intermediate state (where L ^ but 
P — M — 0). They are those which display an antisymmetric linear magnetoelectric 
effect and contain the combined symmetry operation of space inversion followed by 
time-reversal: 

1', 2/m', 2'/m, m'mm, 4/m', 3', 6/W, A/m'mm, 3'm, and 6/m'mm .(4) 

These are all possible space groups of the targeted antiferromagnetic paraelectric 
phase, where a polar distortion can induce a weak magnetization as a secondary order 
parameter. In order to construct the possible polar subgroups we simply add a polar 
distortion along one of the directions connecting antisymmetrical components of the 
linear magnetoelectric tensor to all these groups, and calculate the resulting symmetries. 
This results in the following 5 point groups: 

1, 2', m, m', and 2'm'm . (5) 

We note that these are precisely those determined by Fox and Scott by considering 
all magneto-electric point groups which allow both a spontaneous polarization and a 
spontaneous magnetization and requiring P _L M [13] . 

We point out that even though within our analysis we have first proceeded from 
the paramagnetic-nonpolar case to the non-polar antiferromagnetic case, and then 
subsequently to the polar-magnetic case, it is not required that in the real system 
the magnetic phase transition occurs at a higher temperature than the ferroelectric 
transition. The important point is that an antiferromagnetic paraelectric reference 
phase with the required symmetry can be constructed in principle. In fact displacive 
ferroelectric phase transitions often occur well above room temperature, and therefore 
the critical temperature for the effect described in this article is expected to be 
determined by the antiferromagnetic Neel temperature in most casesO The Neel 

+ In this case, the magnetization will appear at Tjv, rather than at the ferroelectric transition, but the 
effect is still ferroelectrically-induced. 
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temperature of antiferromagnetic oxides with strong superexchange interaction is often 
above room temperature and thus the described effect is not limited to low temperatures. 

6. Summary and Outlook 

In this article we have illustrated the general principle behind ferroelectrically-induced 
ferromagnetism through the DM interaction, focusing on the corresponding symmetry 
requirements. We have discussed in detail why these requirements are fulfilled for R3c 
FeTiOs, where the magnetic Fe 2+ cations occupy the perovskite A sites, but not for 
isostructural BiFeC>3, where the magnetic Fe 3+ cations are located on the perovskite 
B sites. As discussed in Sec. HI the crucial difference is that in the R3c paraelectric 
reference structure inversion centers are located at the midpoints between the A sites 
(structural criterion) but not between the B sites. These inversion centers are destroyed 
by the ferroelectric displacements, and since they are the only symmetry operations 
prohibiting a DM interaction between the A sites (magnetic criterion), the ferroelectric 
distortion "switches on" weak ferromagnetism in FeTiOs. As a result, polarization and 
magnetization in FeTiOs are coupled via the trilinear invariant, Eq. ([3]). 

In Sec. [5] we have reformulated the design criteria of [H] as symmetry requirements 
for the antiferromagnetic order parameter L, and we then constructed magnetic point 
groups that are compatible with ferroelectrically-induced weak ferromagnetism. In the 
most common case, where space inversion is a symmetry element of the paraelectric 
reference structure, this requires L to transform like a magnetic toroidal moment, i.e. 
L has to be odd under both space and time inversion but invariant under the combined 
operation. 

Of course symmetry analysis gives only qualitative information, i.e., it tells us 
whether or not a certain effect is in principle possible for a given symmetry. On the 
other hand, it is also desirable to subsequently quantify the corresponding effect. Here, 
first principles calculations using density functional theory represent an invaluable tool. 
Within certain constraints, these calculations can provide very reliable quantitative 
information about e.g., structural parameters, lattice instabilities, magnetic coupling, 
and magnetization. First principles calculations have been used in [T6| [TTJ [18] to verify 
and quantify weak magnetic order and its relation to structural distortions, both polar 
and non-polar. Combined symmetry analysis and first principles calculations have 
identified FeTi03 (and the related Ni and Mn compounds) as prime candidate for the 
realization of electric-field switchable weak ferromagnetism close to room temperature 
|18j . which now awaits experimental verification. 

Symmetry guided first principles design of novel promising materials is in line with 
the general strategy for a rational computational materials design outlined by Spaldin 
and Pickett [12]: after a candidate material with the right symmetry has been found 
on symmetry grounds, first principles calculations can be used to verify whether this 
material really shows the desired effect and to determine the relevant quantities. 

Such a strategy is particularly useful for the design of novel complex oxide materials, 
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where intrinsic properties can easily be hidden by extrinsic effects, such as oxygen 
deficiency, micro-crystallinity, or the presence of small amounts of competing phases, and 
can thus easily be missed in a purely experimental approach. We therefore believe that a 
rational materials design based on first principles calculations indeed represents a very 
powerful approach to search for novel materials with unexpected and technologically 
useful properties. 
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